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Abstract. We compare our previously proposed hard-thermal-loop (HTL) resummed calculation of quark 
number susceptibilities using a self-consistent two-loop approximation to the quark density with a recent 
calculation of the same quantity at the one-loop level in a variant of HTL-screened perturbation theory. 
Besides pointing out conceptual problems with the latter approach, we show that it severely over-includes 
the leading-order interaction effects while including none of the plasmon term which after all is the reason 
to construct improved resummation schemes. 



1 Introduction 

In view of the ongoing search for quark-gluon plasma sig- 
nals in the early stages of ultrarelativistic heavy-ion colli- 
sions, quark number susceptibilities (QNS) have recently 
received enhanced attention because of their direct con- 
nection with fluctuations of conserved charges which could 
in principle discriminate against a purely hadronic phase 
Concurrently, new results for QNS have become 
available from lattice gauge theory [1|J which consider- 
ably improve upon previous studies [ 3] , and moreover ex- 
tend them to higher T/T c . The diagonaQ QNS are found 
to increase sharply at the deconfinement phase transition 
toward a large percentage of the ideal gas value %o = 
NT 2 /3 for SXJ(N) and massless quarks. 

Conventionally resummed perturbative results for ther- 
modynamic quantities, on the other hand, do not seem to 
be applicable to account for the observed deviation from 
ideal-gas behaviour because of a complete lack of conver- 
gence for all temperatures of interest. In the case of the 
free energy, this problem is particularly severe, because 

3/2 

the so-called plasmon contribution oc a s is larger than 
the leading-order interaction term oc a s for all tempera- 
tures T < 10 5 T C . A similar but less dramatic problem also 
occurs with the perturbative result for the diagonal QNS 



X_ 

Xo 



= 1-2- 



7T 



+ 8^/l + ^(^) 3/2 + 0(^log(a s )) (1) 



1 Offdiagonal QNS are strongly suppressed in the high- 
temperature phase. Still, they give rise to a puzzling discrep- 
ancy between recent analytic and lattice calculations. The 
leading-order effect oc Q 3 log(a s ) has recently been calculated 
in Ref. implying a numerical value of ~ 10~ 4 , whereas the 
available lattice results are claimed to be consistent with zero 
within an accuracy of < 10 -6 m. 



for QCD (N = 3) with N f quark flavours. For N f = 2 
the plasmon term overcompensates the term oc a s for all 
temperatures T < 40T C , and only for T > 700T C does x/xo 
show the expected growth with temperature, starting from 
values extremely close to the ideal-gas result. 

In Ref. M we have shown that these problems can 
be avoided by a reorganization of perturbation theory 
which is based on a self-consistent (^-derivable) two-loop 
approximation to the thermodynamic potential |^]. The 
latter leads to a nonperturbative expression for entropy 
and quark density which can be used to resum the so- 
called hard thermal loops (HTL) and particular 
next-to-leading order corrections thereof. The results for 
QNS thus obtained are monotonic functions of T/T c which 
account at least for a sizeable part of the deviation from 
the ideal-gas behaviour observed in lattice calculations for 
T/T c > 2T C . 

Recently, a different approach to resum the effects of 
HTL in QNS has been put forward in Ref. |ll| which 
starts from quark number charge correlators. Employing 
HTL propagators and vertices at one-loop order, one finds 
substantially larger deviations from the ideal-gas limit, 
seemingly in a good agreement with the lattice results 
ofRefs. j§. 

In view of the large efforts invested at present by lattice 
gauge theorists to explore effects of small chemical poten- 
tials at high temperature in QCD, we think it worthwhile 
to explain the fundamental differences between our ap- 
proach and that of Ref. pl| and why, in our opinion, the 
results of the latter are actually misleading. In particular 
we show that the one-loop results of Ref. severely over- 
include the leading-order interaction effects, while they 
contain none of the plasmon effects oc ct'J 2 (which are 
the source of the problems with conventionally resummed 
perturbation theory). Moreover, we point out a certain 
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technical difficulty that has been overlooked by the au- 
thors of Ref . |ll[] , but has the effect to render their result 
ill-defined in a distributional sense. 

More importantly even, we comment on a conceptual 
problem with the approach followed in Ref. which 
arises because the HTL action is no longer used as the ef- 
fective theory for soft modes, but is used throughout all of 
phase space. Implicitly the definition of the quark number 
charge operator is modified such as to no longer conform 
with the operator employed in lattice calculations. 

Before discussing the approach of Ref. 11 1 in detail in 
Sect. ||, we briefly review the QNS as obtained from HTL- 
resummed thermodynamic potentials. Sect, ^summarizes 
our conclusions. 



2 QNS from resummed thermodynamic 
potentials 

2.1 Generalities 

The QNS of a given quark flavour is by definition the re- 
sponse of the quark number density AT to an infinitesimal 
variation of the associated chemical potential fi, 



dp, 



d 2 P 
dp 2 



/ d 3 x(p(0,x)p(0,0)) (2) 



where P = ((3V)~ 1 logZ is the thermodynamic pressure, 
/3 = T^ 1 and p — ipj°ip. 

When thermodynamic consistency is automatic, for ex- 
ample in strict perturbation theory to a given order in a s , 
it does not matter which of the equivalent expressions on 
the right-hand side of (^) is employed. However, when fur- 
ther resummations are performed that amount to a partial 
inclusion of higher-order effects, it does in fact matter. To 
set the stage we begin by briefly reviewing the approaches 
which focus on the thermodynamic potential before turn- 
ing to Ref. [ p"T| which starts from the quark number charge 
correlator. 

Expressed as a functional of full propagators (D for 
gauge bosons and S for fermions, and assuming a ghost- 
free gauge choice) the thermodynamic potential fl = —PV 
= -T log Z has the form fH 

(3f2[D,S] = -TrlogD-i-^TrllD 

-Tr log S- 1 +Tr£S + $[D, S] (3) 

where <!> is the sum of 2-particle-irreducible "skeleton" di- 
agrams whose lowest-order (2-loop) contributions are 



$[D,S] = -1/12 




-1/8 




+ 1/2 iA/WV\i + . 



As a functional of D and S, fi is subject to the sta- 
tionarity condition, 



8f2[D, S]/SD = = 5Q[D, S]/SS, 



which is equivalent to 

6$[D,S\/8D = ±II, 



5<P[D,S]/6S = S, 



(5) 



for the self-energies 77 and S. Expressing 77 = D 1 — D 
and £ = S~ x — Sq 1 in terms of bare propagators Dq 
and So, the representation (^|) of course reproduces the 
ordinary loop expansion. 

For example, the leading-order interaction terms oc a s 
are given by the 2-loop diagrams in <P, whereas single pow- 
ers of the self-energy insertions in a propagator cancel out 
in the first four terms of the right-hand side of eq. (H) . 

Ordinary perturbation theory, however, has infrared 
problems at finite temperature if the repeated self-energy 
insertions contained in the term ^Tr log D~ x are expanded 
out perturbatively. These can be remedied by a resumma- 
tion of the leading order Debye mass 



m 2 D = {2N + N f ) 



9 2 T 2 



E 



9h± 
2ir 2 



(G) 



in the (chromo-)electrostatic propagator, where g 2 = 47ra s 
(though new infrared problems arise at order a 3 ). Ex- 
panded in powers of g, the resummation of the Debye 
mass in ^Tr \ogD~ 1 gives rise to the so-called plasmon 
term in the pressure 



P 3 = N 3 Tm%/(\2TT). 



(7) 



It is this term which is responsible for the dramatic deteri- 
oration of the apparent convergence of a perturbative ex- 
pansion of P in g at finite temperature, and, as remarked 
in the introduction, to a somewhat lesser degree for QNS 
which can be derived from the pressure. 

2.2 Screened (HTL) perturbation theory 

The loss of apparent convergence upon inclusion of the 
plasmon term in the pressure is in fa ct g eneric and also 
occurs in a simple scalar ip A theory iQ. This problem 
arises as soon as finite-temperature contributions are ex- 
panded out in powers of the coupling, which is necessary 
for the standard ultraviolet renormalization programme 
to become applicable. In order to avoid this, it has been 
proposed uM to reorganize perturbation theory by adding 
screening masses to the classical Lagrangian and to sub- 
tract them as counter-terms, but in contrast to the usual 
resummation programme at finite temperature ^,|TTJ, this 
is done for both hard and soft momentum regimes. This 
in fact alters the ultraviolet structure of the theory, but 
when combined with a simple minimal subtraction of the 
additional divergences this resummation appears to sig- 
nificantly improve the apparent convergence of thermal 
perturbation th eory . 

In Refs. Jll|[l6j this approach has been extended to 
QCD at one-loop level. It amounts to keeping only the 
logarithmic terms in (^|) and replacing D and S by the 
HTL propagators, 



(4) 



p D l-loop-HTL _ Tr log £-l _ Xr lQgS - 



(8) 
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where hatted quantities refer to HTL. If the thermal mass 
parameter in the HTL are exactly the lowest-order ones, 
this includes the correct plasmon term (Q) without causing 
the pressure to exceed the ideal-gas value. However, the 
leading-order interaction pressure cx a s is over-included 
by a factor of 2.[] 

In order to have both, the leading-order term cx a s and 
the plasmon term o^J 2 included correctly, it is necessary 
to go to two-loop order. Starting from two-loop order, one 
can turn this so-called "HTL perturbation theory" into 
a variational perturbation theory, where the HTL action 
is no longer used as an effective theory for soft modes 
as in standard HTL resummation P,|lCj||, but just as a 
gauge invariant mass term which is then eliminated by a 
principle of minimal sensitivity. 

Because of the HTL action involves non-local self-ener- 
gies and vertices, this optimization of perturbation theory 
is extremely difficult and has only recently been carried 
through for QCD to 2- loop order 0. The results are a 
clear improvement over conventionally resummed pertur- 
bation theory as the resummed pressure remains below 
the ideal-gas limit despite a full inclusion of the contribu- 
tions through order a 2 log(a s ). However they appear to 
account for less than half of the deviation from ideal-gas 
behaviour that is observed on the lattice. 

The main virtues of this approach are its systematic 
nature and manifest gauge invariance. From a physical 
point of view, a possible weakness is that the HTL action 
is used uniformly for soft and hard momenta, whereas 
the HTL action is accurate only for soft momenta, and 
for hard ones only in the vicinity of the light-cone. A re- 
lated problem is that the artificial UV divergences that are 
introduced involve new subtraction scheme dependences. 
While these start to be suppressed by powers of a s only 
at the (rather forbidding) three-loop order Q, these ad- 
ditional scheme dependences turn out to be numerically 
rather weak in the two-loop result for QCD. 



2.3 HTL-resummation of the 2-loop ^-derivable 
entropy and density 



While HTL-screened perturbation is in principle rather 
generally applicable, we have found, following up an ob- 
servation made in Ref. |19"| , that specifically for the first 
derivatives of the thermodynamic potential one can de- 
rive remarkably simple expressions from a self-consistent 
2-loop approximation to the skeleton expansion (||) of the 
QCD thermodynamic potential. Because of the station- 
arity property, these derivatives act only on the explicit 
statistical distribution functions, and not also on those 
contained in propagators and self-energies. Moreover, af- 
ter differentiation, the contribution from <p 2 - loo P just can- 
cels part of the second and fourth term on the right-hand- 
side of (0). The derivatives with respect to temperature 



Fig. 1. Next-to- leading order corrections to the asymptotic 
fermion mass 



and chemical potential give entropy and quark densities, 
respectively, reading (^] 

5 = / S^r^ [Im log.0- 1 - Im 77 Re D] 



N = 



-2tr 
-2tr 



(2tt) 4 dT 
d i k df 

(2tt) 4 df 
d 4 k df 

(2tt) 4 c^ 



Im logS- 1 - Im^ReS , (9) 
[Im logS 1 " 1 - JmE Re S] , (10) 



2 As explained in the last paper of Ref. ||, Ref. [|L5| had 
an even stronger over-inclusion due to an inconsistent use of 
dimensional regularization. 



where D and S are determined by one-loop gap equations 
(||) obtained by restricting <P to two-loop order. 

Although these equations have the form of one-loop 
expressions involving dressed propagators, they include 
all the two-loop contributions, but incorporated in the 
spectral properties of the quasi-particles described by the 
dressed propagators. This implies that both the leading- 
order interaction terms cx a s and the plasmon effect cx 
aj" 1 are completely taken into account as soon as D and 
S are evaluated to sufficient accuracy. 

Because the expressions (^|) and ( |To| ) are manifestly 
ultraviolet finite, they can be used to resum the effects of 
HTL without the necessity of subsequent expansions and 
truncations.^ At soft momenta, the HTL are valid expres- 
sions to the actual full propagators that one would have to 
use in a self-consistent scheme; at hard momenta it turns 
out that to leading order the above expressions only probe 
the vicinity of the light-cone where HTL self-energies re- 
main accurate. The next-to-leading order effect, which is 
the plasmon effect, turns out to be to one part covered 
by HTL resummation in the soft regime, and to the re- 
maining part by corrections to the so-called asymptotic 
thermal masses from HTL-resummed one-loop diagrams. 

In the case of the quark number density functional, 
from which the QNS can be derived, it turns out that all 
of the plasmon effect is associated with next-to-leading 
order corrections to the asymptotic fermion mass, shown 
in Fig. 

If only the HTL approximation to the fermion propa- 
gator is employed, the quark number density still contains 
the complete leading-order interaction effects cx a s , and 
- since it need not be expanded out perturbatively - also 
subsets of higher-order effects, namely those associated 
with repeated HTL insertions. 

In Ref. [7| we have evaluated the QNS obtained by 
taking the derivative of ( |l0| ) with respect to /i, both in 
the HTL approximation and in a next-to-leading approxi- 
mation which incorporates the plasmon effect through the 
corrections to the asymptotic fermion mass from the dia- 
gram of Fig. [I]. 

3 For a different approach based on the pressure see also 
Ref. H. 
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Fig. 2. Lowest-order contributions to (|l^) in bare perturbation 
theory. 



3 QNS from HTL-resummed charge 
correlators 

In Ref. |ll| an HTL-resummed QNS has been constructed 
by starting from the charge correlator 

X = P f d 3 x(p(0,x)p(0,O)) = B [ d 3 z77>(0,x) (11) 



where IV? V is the current-current correlator of a given 
flavour charge (suitable linear combinations of such quan- 
tities give the correlators of electric charge and baryon 
number) . 

In Fourier-space one has 21 , 22], p3]] 



X 



lim 3 

fc->o 



Zir 



(12) 



where the well-known fluctuation-dissipation theorem al- 
lows one to write (in the notation of Jul) 



(13) 



ujn> 



with n^ v the retarded response function. 

In conventional perturbation theory the first few dia- 
grams contained in (jl^) are shown in Fig. The first 
diagram gives the ideal gas value, and the leading-order 
interaction term oc a s in ([l]) is given by the two-loop or- 
der diagrams. The plasmon effect oc o?J 2 comes from those 
higher-order diagrams which correspond to repeated self- 
energy insertions into the gluon lines of the two-loop dia- 
grams of Fig. 0. 

Calculating \ from (Q) is in fact a bit more involved 
than starting from the thermodynamic potential as also 
noticed in Ref. Because charge conservation implies 
(u>, 0) = one has linifc^o IIq (uj, k) ~ S(u>). 
If, for example following HTL-screened perturbation 
theory, one is interested in the one-loop contribution aris- 
ing from dressing the propagators in Fig. ||, this will spoil 
this behaviour. In order to have charge conservation one 
needs to employ HTL vertices in addition to HTL prop- 
agators, which is precisely what the authors of Ref. |yj 
have proposed. 

However, this raises an important conceptional prob- 
lem: in effect this use of the HTL vertices replaces the 
ordinary charge operator tpj°ip by the non-local object 
■^(j +r°)ip derivable from the non-local HTL action. The 
correspondingly re-defined QNS is therefore no longer di- 
rectly related to the quantity defined in ( |lT| ) and measured 
in lattice simulations. 



4 One-particle-reducible diagrams, which in principle con- 
tribute, vanish because of the tracelessness of colour matrices. 




+2- 






Fig. 3. One-loop contributions to ( |l^ ) in HTL-screened per- 
turbation theory. The vertex parts built from bare propagators 
are understood to be evaluated in the HTL approximation. 



-o- 





+2 



+ <E> 



Fig. 4. Two-loop contributions to (|12|) in HTL-screened per- 
turbation theory which contribute to the leading-order inter- 
action terms oc a s . 



This is in fact a problem only because the HTL action 
is no longer used as an effective theory appropriate for soft 
momentum scales, but is used equally for soft and hard 
momenta. As an effective theory, obtained after integrat- 
ing out the hard momenta and used for soft modes only, 
the appropriate charge operator is indeed the non-local 
quantity involving the HTL vertex r°. It is this opera- 
tor which enters in a perturbative matching to the full 
theory. But replacing the ordinary charge operator by the 
HTL-dressed one for all momenta clearly corresponds to 
abandoning the definition ( |ll"|) for the QNS. 

Leaving this issue aside for now, we continue by analy- 
sing the diagrammatic content of x m HTL-screened per- 
turbation theory. Since HTL vertices are strictly one-loop 
quantities, to one-loop order it is as shown in Fig. ^. 

However, while all the topologies that are present in 
the two-loop diagrams of Fig. also appear in the dia- 
grams of Fig. ||, their combinatorial factors are different. 

This shows that in the one-loop HTL approximation 
to ([I2]) the a s -contributions are all overcounted. The sec- 
ond diagram of Fig. || is contained with correct combina- 
torics in the first diagram of the right-hand side of Fig. |^ 
but appears another time through the HTL 4-vertex; the 
third diagram of Fig. |^ is seen to be over-included by a 
factor of 2. Moreover, because the HTL approximation 
for the undressed self-energies and vertex subdiagrams in 
Fig. |^ does not provide the complete leading-order terms 
for hard inflowing momenta, this leads to a further source 
of incompleteness of the terms oc a s . 

The authors of Ref. jll| do not specify how to ex- 
tend their approach to two-loop order. It is clear, how- 
ever, that the correct counting is restored only when the 
modification of the quark number charge is undone by 
HTL counter-terms and all two-loop diagrams are added. 
The relevant diagrams for completing the order a s result 
are shown in Fig. W, where the first and the third diagram 
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Fig. 5. x/xo as a function of M/T. The dashed line is the 
perturbative result to order g 2 , the upper full line is the 2-loop 
^-derivable approximation evaluated with HTL propagators, 
the lower full line is the 1-loop HTL result of Ref. 111. 



correspond to HTL counter-terms to the charge operator .|j 
This means that the definition of the charge operator in 
( |ll| ) has to be modified order by order to approach to 
standard definition of QNS at least at infinite loop order. 

3 /2 

This also shows that the plasmon term oc a s , which 
is the reason for seeking improvements of conventionally 
resummed perturbation theory, only appears in the two- 
loop order diagrams of HTL-screened perturbation theory, 
namely through the dressed vector boson lines with a blob. 
The vector boson lines within the HTL vertices of Fig. || 
are not dressed and thus do not capture anything of the 
plasmon effect. 

We are now in a position to compare with the HTL- 
resummation approaches discussed in the previous section. 

In one-loop HTL-screened perturbation theory along 
the lines of Refs. |T^,[l6j the leading-order interaction term 
to the thermodynamic potential is over-included by a fac- 
tor 2, but the plasmon effect is complete (as long as only 
the leading-order HTL mass parameter is used; including 
higher-order corrections in the latter would spoil this). 
The QNS have not been calculated in this approach, but 
the same pattern would apply, provided the perturbative 
HTL masses are inserted before differentiating with re- 
spect to 

On the other hand, in the 2-loop ^-derivable quark 
density (|lfj|), the leading-order interaction term is already 
correctly included when evaluating it in the HTL approx- 
imation but the plasmon term is absent (it arises exclu- 
sively from next-to-leading order corrections to the asymp- 
totic (hard) thermal fermion mass). 

In Fig. |, the QNS we have obtained from ( |To| ) in the 
HTL approximation [Q is evaluated as a function of the 



5 The 2nd and 4th diagram have opposite combinatorial fac- 
tors but may also contribute to order a s because of the incom- 
pleteness of the HTL approximation for hard loop momenta. 

6 Giving up thermodynamic consistency, one could think of 
identifying the mass parameter in HTL-screened perturbation 
theory only after this differentiation. This would in fact give 
a correct leading-order interaction term, but would loose the 
plasmon effect (and not completely reproduce the terms of or- 
der a 2 and higher contained in (|1C|)). 
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Fig. 6. R = (x - Xo)/{x; ' th . - in) as a function of M/T 
for the 1-loop HTL result of Ref. (ll) (upper half of the plot) 
and for the 2-loop ^-derivable approximation evaluated with 
HTL propagators (lower half of the plot - note the different 
scale!). The perturbative result to order g 2 corresponds to the 
value 1. Only the 2-loop HTL result approaches 1 in the limit 
M /T — > 0, where eventually perturbation theory should be 
reproduced; the 1-loop HTL result of Ref. |ll[] is seen to over- 
include the leading-order interaction effect by a factor which 
diverges logarithmically as M/T oc g — > 0. 



fermionic plasmon massQ M/T and compared with a nu- 
merical evaluation of the one-loop HTL result reported in 
Ref. JO]. While our result shows slightly slower deviation 
from the ideal-gas limit than the strictly perturbative re- 
sult to order a*, the result of Ref. |llj has considerably 
stronger deviations because of the over-inclusion of the 
leading-order interaction term. 

In Fig. [| we consider the expression 



R=(X~ Xo)/(Xp.l h . - Xo) 



(14) 



which measures the deviation of the interaction part of \ 

(2) 

from the perturbative result X p th ^° or der a s , and plot 
the respective results, again as a function of M/T. While 
our result obtained from ( |l0|) in the HTL approximation 
goes to 1 in the limit of a weakly coupled theory, the ef- 
fective over-inclusion of the leading-order interaction term 
of Ref. pd[ | diverges in this limit. In fact, it turns out that 
the result reported in Ref. [pd| involves a contribution 
oc (M/T) 2 \og(M/T) ~ a s \og(a s ), which docs not exist 
in the correct perturbative expansion. This over-inclusion 
problem is therefore much more severe than in the case 
where 1-loop HTL-screened perturbation theory is applied 
to the thermodynamic potential 15,19]. 



For plots in terms of a 3 or T/T c see Ref. 
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This clearly shows that the one-loop HTL resumma- 
tion of the charge-charge correlator cannot be compared 
with either our results, which are based on the two-loop 
expression ([n]), or perturbation theory which it seeks to 
improve upon. In our opinion it is also completely prema- 
ture to compare with the available lattice results on QNS, 
because the two-loop contributions of Fig. |I] will have to 
correct for the enormous over-inclusion of terms oc a s . 

But it appears to be questionable whether a two-loop 
HTL-screened perturbation theory calculation of ( [l2| ) is 
at all practicable. There are in fact certain technical prob- 
lems with the result reported in Ref. already at one- 
loop order. In eq. (34) of Jll| one can see that the result 
for IIq (uj, 0) is proportional to the integral 



d 3 k J dx J dx'nF(x)n,F(x')p-i r (x 1 k)p^(x,k) 
(uj — x — x') 2 S(lu — x — x') 



(15) 



where p±{x, k) are the HTL spectral functions for the two 
fcrmionic quasi-particle branches of the HTL approxima- 
tion. In Ref. |[l| the latter are used to put x — —x', so that 
the second line of ( |l5| ) is reduced to S(lo) in conformity 
with the expectations from charge conservation. However, 
this term is clearly ill-defined and might with equal justi- 
fication be put to zero identically as is suggested by the 
way we have written it. In order to have a well-defined 
expression, it seems to be necessary to keep the external 
spatial momentum different from zero and take the limit 
to zero only after having performed the integral over w (as 
demanded by (|l2[)). But that would make its evaluation 
in HTL perturbation theory a hopelessly difficult task, al- 
ready at one-loop order. 



4 Conclusions 

In this paper, we have discussed various possibilities for 
HTL-resummation in the calculation of QNS and have in 
particular analysed the recent proposal of Ref. |p!l| . We 
have shown that the resummed one-loop calculation pre- 
sented there severely over-includes the leading-order inter- 
action terms, while not including anything of the plasmon 
effect, both of which would be corrected only at two-loop 
order. Thus only the latter should be viewed as an im- 
provement over ordinary perturbation theory and used as 
such in a comparison with the available lattice results. 

However, because of the technical problems mentioned 
at the end of the previous section, it would seem to be 
more sensible to calculate the QNS through a 2-loop HTL- 
screened perturbation theory evaluation of the pressure 
along the lines of Ref. [jl7| . 

On the other hand, the HTL-resummed calculation of 
QNS of Ref. § is based on a 2-loop ^-derivable approxi- 
mation and does include correctly both the leading-order 
interaction effect, ~ a Sl and the next-to-leading-order one, 
~ a'J 2 (together with an infinite series of higher order ef- 
fects due to HTL). The results in Ref. flTtl show the same 



trend as the lattice results, but a significant difference still 
remains, which calls for further studies, both on the ana- 
lytic side, by further improving the resummation schemes, 
and on the lattice side, by increasing the reliability of the 
numerical results. 
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